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Abstract 

A new method of refrigeration is proposed. Cooling is obtained by therm- 
ionic emission of electrons over periodic barriers in a multilayer geometry. 
These could be either Schottky barriers between metals and semiconduc- 
tors or else barriers in a semiconductor superlattice. The same device is an 
efficient power generator. A complete theory is provided. 
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A new method is proposed for refrigeration and power generation. The 
devices are composed of multilayers of periodic barriers. The currents are 
perpendicular to the barriers. Such devices are presently under consideration 
for thermoelectric cooling and power generation. Here we propose that the 
same devices can be used for efficient thermionic cooling, or power generation. 

What is the difference between a thermoelectric device Jl], [| and a therm- 
ionic one|| 10]? Their descriptions are remarkably similar. In both cases 



one puts a temperature gradient on a semiconductor. Both devices are based 
upon the idea that electron motion is electricity. But the electron motion 
also carries energy. Forcing a current transports energy for both thermionic 
and thermoelectric device. The basic difference seems to be whether the 
current flow is ballistic or diffusive. In thermionic motion, the device has 
relatively high efficiency if the electrons ballistically go over and across the 
barrier. They carry all of their kinetic energy from one electrode to the 
other. In thermoelectric devices, the motion of electrons is quasiequilibrium 
and diffusive. One can describe the energy transport by a Seebeck coefficient, 
which is an equilibrium parameter. There have also been recent theories of 
nonequilibrium thermoelectric effect s [pJj1 . They seem to be midway between 
the thermionic(ballistic) and thermoelectic (quasiequilibrium) regimes. We 
note that many semiconductor superlattices, with short periods, are being 
made [[12] and measured along the c-axis. These periods are so short that 
electron motion over the barrier is probably ballistic. We suggest that these 
devices cool by thermionic effects rather than thermoelectric ones. 

Earlier we suggested that thermionic devices could be used as refrig- 
erators. For the usual device of two metal plates separated by an air gap, 
cooling at room temperature requires a low work function. No metals have 
values that low (0 ~ 0.3 eV). Numerous groups have suggested [3 -jTOj that 
such small barriers are easily attainable in semiconductor systems. Here the 
barriers are semiconductors, while the electrodes could be either metals or 
other semiconductors. We considered this geometry in our original paper, 
but thought that the thermal conductivity of the solid would be a major 
obstacle. Here we show that the thermal effects can be dealt with by going 
to a multilayer geometry. We show that such devices have efficiencies which 
could be twice those of thermoelectric ones. 

The physics behind thermionic cooling is simple. Most physicists are 
familiar with the technique of cooling liquid helium-4 by pumping the vapor 
from the cryostat. The most energetic helium atoms leave the liquid and 
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become gas molecules. Pumping them away removes these energetic atoms, 
thereby cooling the liquid. In thermionic refrigeration, one uses a voltage 
to sweep away the most energetic electrons from the surface of a conductor. 
Those electrons with sufficient energy to overcome the work function are 
taken away to the hot side of the junction. Removing the energetic electrons 
from the cold side cools it. Charge neutrality is maintained at the cold side 
by adding electrons adiabatically through an ohmic contact. 

Thermionic devices must have the electrons ballistically traverse the bar- 
rier in order to have a high efficiency. This requires that the mean-free-path 
A of the electron in the barrier be longer than the width L of the barrier. 
This constrains the barrier width L to be rather small. This fact is a key 
feature of the analysis. The general constraints are 



where L t is the minimum thickness to prevent the electron from tunneling 
through the barrier. Values of L t seem to be around 5-10 nm for most semi- 
conductors. The values of A are less certain. There are good measurements 
of the electron mean-free-path for motion along the layers, but little data for 
motion perpendicular to the layers. If these values are similar, then one can 
easily find that A ~ 50-100 nm. Thus there is room in the above inequality 
to find values of L which work. 

The general form of the energy currents for a single barrier are 



The first term in the energy current is the electron part from thermionic emis- 
sion. It is given below. The second term is the phonon part, which contains 
the thermal resistance R\ for one barrier. The thermal resistance depends 
upon the thickness (L e ) and thermal conductivity (K e ) of the electrodes, as 
well as two interface ('Kapitza') terms -R/|13|, 0]. The largest term will usu- 
ally be L/K for the semiconductor barrier. The problem is that L is small 
which makes R\ small which makes 5T/Ri big. This is the problem with the 



A > L > L, 



(1) 




(2) 




(3) 



(4) 
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thermal conductivity. It can be overcome by having 5T be small. This is the 
reason for the multilayer geometry. Each barrier can support only a small 
temperature difference 8%. A macroscopic temperature AT is obtained by 
having N layers so that AT = N(STi). Most of our modeling assumes that 
5Ti has an average value of 1-2 °C. Below we show that adequate cooling 
power is available for these values. 



1 Refrigeration 
1.1 Single Barrier 

The first step in the derivation is to solve for the currents over a single barrier 
of width L. We assume the barrier is a constant at zero applied voltage, 
which means it has a square shape. At a nonzero voltage 5V the barrier has 
the shape of a trapezoid. The formulas for the electrical (J) and heat (Jq) 
currents are given in terms of the hot and cold temperatures (T^TjjS] 

(5) 

(6) 
(7) 

le<t> + 2k B T h ]J Rh e-* sv ^- 6 -^- (8) 

Hi 

ST = T h -T c (9) 

These equations are for a single barrier. Eqn.fl^) is the standard Richardson's 
equation[15|] for the thermionic current over a work function ecf) which in this 



Jrj = 


AT 2 e - e <t>/ k BTj 


A = 


emk% 




2n 2 h 3 


J = 


Jrc — JrhZ 6 


eJq = 


[e(p + 2k B T c ]J Rc 



case is the Schottky barrier height between the metal and semiconductor. 
Alternately, it is the barrier in a semiconductor quantum well. The factor 
of T denotes the fraction of electrons transmitted from the metal to the 
semiconductor. It is calculated using quantum mechanical matching of the 
wave functions. The formulas for J and Jq assume the bias eSV is to lower 
the Fermi level on the hot side, so that the net flow of electrons is from cold 
to hot. 

For electrons the charge e is negative which makes 0, J also negative. We 
find this confusing to treat, so we take e and <fi as positive as if the system 
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were a hole conductor. Then J > for particle flow to the right, which we 
are assuming. 

We show that for a single layer the optimal value of applied bias eSV oc ST 
which is also small. Denote as T the mean temperature of the layer, and then 
T c = T — ST/2, Th = T + ST/2. Then we expand the above formulas for the 
currents in the small quantities (ST/T, eSV/kBT) and find, after some algebra 

J = ^IW-Vj] (10) 
J Q = J R (b + 2)[SV -V Q ] (11) 

" - £ « 12 » 

eVj = k B ST[b+2] (13) 
eV Q = k B STq (14) 
q = b + 2 + u (15) 

u = TT~o" ( 16 ) 



6 + 2 

e 

k B R\JR 

ek B (T R 



Z e b (17) 

(18) 



RiA(k B T) 2 T V T 
(ksT R f = -^L- (19) 

The symbol J/? denotes Richardson's current eqn.(|^) at temperature T. The 
dimensionless constant Z plays an important role in the results. Refrig- 
eration requires that SV > Vq. High efficiency requires that Zq < 1. We 
have also introduced the dimensional constant Tr which is determined by 
the thermal resistance and the effective mass. 

The efficiency r\ of a single barrier is defined as the heat flow divided by 
the power input 

- - Jq (20) 



SVJ 

k B T{b + 2) 5V- Vq 
e SV(SV - Vj) 



(21) 
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We now vary 5V in the above equation to find the value of 5V = V m which 
gives the maximum efficiency r] m which is 



V m = Vq + y Vq(Vq — Vj) (22) 

w)* (23) 
m = Jttw (24) 

The first factor on the right in eqn. (^) is approximately the Carnot efficiency 
of the one-layer system. The remaining factor fj is the reduction of this 
efficiency due to the inefficiency of the device. The latter factor depends 
upon the parameter b = e^/ksT and the dimensionless factor Z$ contains 
the thermal resistivity R\. Fig. 1 shows the function fj(4>) for several values 
of Tr in the range between 100 and 500 K. It is important to have fj be as 
large as possible. Clearly this constraint limits the value of Tr to smaller 
values, which has Z to be less than one. For larger values of Z the effective 
efficiency becomes too small to be useful. 

The value of Tr is reduced by the Kapitza resistance Rj at the interfaces. 
Numerical estimates for Ri are given by |T3| , [14]] . The device will only be 
efficient when Rj large. Large thermal boundary resistance can be achieved 
by having acoustic mismatch between the metal and semiconductor layers. 
The important point is that the thermal conductivity term has to stay small. 



The thermal resistivity is large in semiconductor superlatt ices |16|, [17], [L8|, [L9 

HP- 



The results for the single barrier demonstrate that one has to have small 
values of barrier height e<p ~ 2/c#T and large values of the thermal resistance. 
These characteristics carry over to the theory of the multilayer device. Small 



values of the barrier heights are reported in several references [O, [2^. |23l El 



1.2 Multiple Barriers 

Now consider the multilayer thermionic refrigerator. There are N barriers, 
with alternating electrodes. Thus there are 2jV+1 layers, assuming electrodes 
come first and last. Denote 5Ti and 5Vi as the temperature and voltage 
change across one barrier. These values change from the initial to the final 
layer. This change is due to the generation and flow of heat. At each barrier, 
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the electron ballistically crosses the barrier region, and then loses an amount 
of energy eSVi in the electrode. The heat generated at each electrode must 
flow out of the sample according to the equation 

= J f < 25 > 

where Lj is the effective width one one barrier plus one metal electrode. We 
take Li to be a constant, although it could vary with i. The width of the 
device is D = NL{. In normal usage there would be a negative sign on the 
right-hand side of this equation. It is absent since we changed the sign of J. 

The functions change slowly with % and we treat them as continuous 
variables 6Ti/Li = dT/dx and SVi/Li = dV/dx. The variation of 5Vi upon 
i is unknown, but we guess that it is proportional to <5Tj. Introduce the 
dimensionless function v(x) 

e5Vi = k B 5T t v (26) 
J = J^v-b-2] (27) 
eJ Q = J R (b + 2)k B 5Ti[v - q] (28) 

We use eqn. (^7]) to define STi since the current J is the only constant among 
all of the parameters 

(29) 
(30) 
(31) 

(32) 
(33) 

The variable I(x) is introduced in eqn.( p9l|30"D . After trying various alterna- 
tives, we decided it is the convenient variable for the present problem. The 
parameter y is a dimensionless current. Eventually we vary this parameter 
to find the optimal efficiency The thermal resistance of the entire device is 



5T = 








uJ R 




u 


I = 






v-b-2 


D d ^ = 


ITyZ 


dx 


u 


y = 


'(t)«* 




Rn = 


NR X . 
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denoted as Rn- Its value is irrelevant since it ocurs in y which is a variational 
parameter. Use eqn. fl29D to evaluate 8Vi in (l26f) and this result is inserted 
into eqn.([Z5|) 

Jq = J(0 + 2^)(l-/) (34) 

D±{I) = -y^[ u + I(b + 2I)} (35) 

dT yITZ 
D— = - 36 

ax u 

Eqns.(|35|) and ([36]) describe the variation in temperature and voltage across 
the multilayer device. The input parameters are the barrier height e<p(b = 
e(f)/kBT) and the parameter Tr which depends on the thermal resistance. 
The parameter y is variational. One solves the equations for different values 
of y and chooses the value which give the best results. One selects an initial 
guess for I(x — 0) = I c at the cold end along with T(0) = T c . Then one can 
iterate the above two equations using one sided derivatives. The value of I c 
is varied until one finds the desired value of 7\. The calculation is done for 
different values of y until the maximum efficiency is attained. The Carnot 
efficiency is defined as the heat taken from the cold side Jq c divided by the 
input power JAV 

Jqc = j(—){b c + 2)[l-I e ] (38) 




AV = E^=(^)E^^[^ + ^ + 2)] (39) 

where AV is the net voltage drop across the multilayer device. The above 
integral for AV is easy after one has solved for T(x) and I(x). 

Fig. 2 shows the efficiency of a refrigerator, with T c =260 K and T^= 300 
K, for four values of Tr = 100 — 500i^. They are shown as a function 
of of 0. The efficiencies are quite high: i]=2.5 for Tr = 100 K to 77=0.6 
for Tr = 500 K. The comparable value for the thermoelectric refrigerator, 
with a dimensionless figure of merit equal to unity, is r] TE = 0.70. That is 
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comparable to the thermionic results when Tr = 500 K. Thus one wants to 
operate the thermionic device with values of Tr less than 500 K. 

The factor of Tr is determined by the thermal resistivity. Reducing the 
thermal resistivity by a factor of 25 only reduces the efficiency by a factor 
of 5. Thus the efficiency scales with Tr rather than with R-y. This is quite 
different than in thermoelectric devices, where the result is almost directly 
proportional to the thermal resistivity. Note in eqn. (|35|) that Z enters in the 
formula in the combination of Zj (2 + Z) which saturates at large values of 
Z. The increase is Zq is offset by a lower value of barrier height <fi. 

The exact results show that the thermionic refrigerator has about 30% 
of the ideal Carnot efficiency AT/T C for values of Tr = 200-300 K. Our 
calculations show that the efficiency slowly declines with increasing values of 
Tr. The value of where the maximum efficiency occurs also declines with 
increasing value of T R . 

2 Multilayer Thermionic Power Generation 

Power generation can be understood by considering a simple device with one 
barrier between two electrodes. If the electrodes are at different tempera- 
tures, and if there is no initial voltage between the electrodes, then electrons 
are thermally excited over the barrier. A net electron flow goes from hot 
to cold. If the electrodes are insulated, then the electrodes become charged 
and the system develops an open circuit voltage, which opposes the flow of 
electrons and reduces it to zero. Connecting the electrodes to an external 
circuit causes current to flow. Power can be extracted from the device. The 
behavior is identical in concept to a solar cell. Here the efficiency is calcu- 
lated using the usual definition: it is the external power JAV divided by the 
heat extracted from the hot electrode Jq/j. 

2.1 One Barrier 

Again the barriers are assumed to be thinner than the mean-free-path of the 
electrons, so that one can apply the formulas of thermionic emission. There 
is a small temperature drop 5T and voltage drop 8V across the barrier. The 
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efficiency of a single barrier using eqns. (|H)|)-(|T4|) is 

, = 6 m5V ~ V -ll (40) 

' k„T(b + 2) SV-Vq y ' 

The voltage 5V is varied to find the value 5V m which gives the maximum 
efficiency r\ m 

sv m = ^^Vv(Vv-V^) (4i) 

Vm = jrV (42) 

where fj is defined in eqn.(^). It is the same function which is found in the 
efficiency of the refrigerator. 



2.2 N-Barriers 

Again assume that the device has 2N + 1 layers. Barriers alternate with 
electrodes, with electrodes on both ends. The electrodes can be metals, or 
else conducting semiconductors. The barriers are going to be nonconducting 
semiconductors. Again we shall find that the barrier height e<j) is going to be 
small. There is a small temperature drop and voltage drop SVi across each 
barrier i. Eqns.(|l0|)-(|l4l) also apply to the generator. The difference from 
the refrigerator is that now J < and Jq < 0. This makes the parameter 
y be negative. Another change is that b + 2 > v so that the function I(x) 
is negative. Thus we can apply eqns.(^,^) with the understanding that y 
and I(x) are both negative. 

The formula for the efficiency is the inverse of the one for the refrigerator 



(e<f> + 2k B T h )(l-I h ) v ' 

These equations were solved on the computer. We assumed that T c = 300 
K and Th = 400 K, so that AT = 100 K. Figure 3 shows the efficiency as a 
function of barrier height for several values of Tr. This should be compared 
with a thermoelectric generator which has an efficiency of r/=0.048 for the 
same operating temperatures. Clearly the thermionic device is more efficient 
for small values of Tr. Figure 4 shows the temperature and 5Vi profiles along 
the thermionic device, assuming that iV = 100. 
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3 Thermoelectric Analogy 



Thermionic devices are not thermoelectric devices. However, it is useful to 
examine the analogy between the two kinds of devices. Once we linearize the 
eqns. (|IT||Tl~| , p5l) for small voltage drops, and then make them continuous in 
a multilayer geometry, they become identical in form to the equations of a 
thermoelectric. This analogy gives the effective formulas for the conductivity 
a, the Seebeck coefficient S, and the thermal conductivity K. 



a = 


eJ R Li 
k B T 


(44) 


s = 


T< 6 + 2 > 


(45) 


K = 




(46) 


Z = 


aS 2 T b + 2 
K u 


(47) 



The first term in the thermal conductivity is the electronic contribution K e 
Note that it obeys a modified Wiedemann-Franz Law 




K e = crT — (48) 



The last line in eqn. (|4"7D gives the dimensionless figure of merit Z. It is 
known in thermoelectric devices that the maximum efficiency is obtained 
by having Z be as large as possible. Varying b to maximize Z gives the 
relationship 

b £ 2mk B R 1 

where T is the mean temperature in the device. This is an accurate estimate 
of the optimal barrier height 0, where b = etfi/ksT. Furthermore, the ther- 
moelectric estimates of the efficiency of the refrigerator (j] r ) and generator 

7T C — Th 

Vr = AT( 7 + 1) (50) 
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V. = (51) 
7 = VTTZ (52) 

are remarkably accurate compared to the computer solutions. The Eqs. ( j50|) 
and fl5"T|) estimate the correct efficiency with an accuracy of about 1%. These 
formulas provide analytical estimates of the efficiency. One can also use the 
thermoelectric analogy to find the current densities at maximum efficiency 
as well as other parameters. 

In thermoelectric devices it is highly desirable but rather rare to have 
values of Z > 1. However, in modeling thermionic devices we find values of 
Z much larger than one for reasonable values of thermal resistance. The ef- 
fective Seebeck coefficient in eqn.fliSD is about 250-300 /xV/K since values of 
b are between one and two. By using the low values of thermal conductivity 
reported along the c-axis of a superlattice, we estimate that Tr = 200-400 K. 
For this range of parameters the effective dimensionless figure of merit Z is 
between 2 and 5. The efficiencies of the thermionic devices are correspond- 
ingly much higher than in thermoelectric devices. Ballistic transport carries 
more heat than diffusive flow. 



4 Discussion 

A detailed theory is presented of the properties of a multilayer thermionic 
refrigerator and power generator. It is showns that if the thermal resistance 
is high, that the devices can be twice the efficiency of the equivalent ther- 
moelectric devices. The values of thermal resistiviity needed to make them 
work well are in the range of reported values for multiple quantum wells. It 
appears these devices are practical. 

One of the interesting questions is to select the value of N. This deter- 
mines the number of multilayers. As long as this number is larger than about 
ten in refrigerators, it does not change the efficiency. The choice of N may 
affect the thermal resistance. For a fixed temperature drop AT, the value 
of iV affects the average temperature drop per layer (5T) = AT/N. The 
cooling power is the energy current from the cold side. An estimate is 

AT(f)AT /t 
Jq c = r; — exp(-e0/ k B T) (53) 
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where T is the average temperature of the device. For a refrigerator take A 
= 120 A/(cm 2 K 2 ) along with f = 280 K, <p = 0.050 eV, and (5T) = IK. 
Then the result is Jq c = 212 W/cm 2 . This is a large cooling power. Thus 
having a small temperature drop permits ample cooling. It does require a 
small barrier height. Incidently, the numerical solutions allow an accurate 
calculation of Jq c and the above formula is an underestimate by about a 
factor of two or three. One could increase iV by a factor of ten and still 
have ample cooling power. The manufacturing costs would depend on N. 
Another issue is mechanical strength. A thicker device (larger N) is stronger. 
Thus there are trade offs between cooling power, manufacturing costs, and 
mechanical strength. Such engineering issues are beyond the scope of the 
present discussion. 

We thank L. Woods for helpful discussions. Research support is acknowl- 
edged from the University of Tennessee, and from Oak Ridge National Lab- 
oratory managed by Lockheed Martin Energy Research Corp. for the U.S. 
Department of Energy under contract DE-AC05-96OR22464. 

Appendix 

Here we examine the transport of electricity and heat using the drift-diffusion 
eqn.(l). We will show that the amount of heat is negligibly small when 
one uses this equation. Unless the layer thickness is in the regime where 
thermionic emission is valid, the semiconductor barrier does not transport 
significant amounts of heat. 

We consider a square barrier for the semiconductor. Assume there is a 
small applied potential and a small temperature difference 5T = T h — T c . 
Then the formula for the current should be 



J 



= ejj,n— jj,k B T 



n(L) -n(0) 
L 



(54) 



ra(0) 
n(L) 




(55) 
(56) 
(57) 



n 



= n e- e ^ kBT 



We expand the density exponents using T^ c — T ± 5T/2 and find 



J = e -^[sv-v D ] 



(58) 
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V D = bk B 5T (59) 



The above formula resembles eqn.(10). The two constant voltages Vj and 



Vd are similar. The prefactors are very different. Denote by r the ratio of 
the two prefactors, which gives the ratio of the magnitudes of the currents 
predicted by drift-diffusion compared to the currents predicted by thermionic 
emission 

r = TJ^F) m 

- - ^ m 

where a is the conductivity of the semiconductor and L is a characteristic 
length which we take to be one micron. At room temperature we find that 
o"o = 8-1 kS/m. InSb has a mobility of 8 m 2 /(V s) at n — 10 20 /m 3 for 
a = 130 S/m. Even for b = 4 this predicts a ratio of r w 1CT 3 . So the 
current due to drift-diffusion is one thousand times smaller than the current 
due to thermionic emission. The energy currents have the same ratio of 
proportionality. Thus the flow of heat is negligible when the thickness of the 
semiconductor is greater than the mean free path of the electron, so that the 
particles diffuse. 
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Figure Captions 



1. The reduction in the Carnot efficiency fj for a single barrier as a function 
of barrier height for values of T R = 100, 200, 300, 400, and 500 K. 

2. Efficiency r] of a multilayer thermionic refrigerator with T c = 260 K 
and T h = 300 K. Results are plotted as a function of barrier height for 
different values of T R . 

3. Exact efficiency of a multilayer thermionic power generator as a func- 
tion of barrier height <fi for several values of T R . We assumed T c = 300 
K and T h = 400 K. 

4. Change in the temperature Tj and voltage Vj as a function of position 
along the multilayer power generator. We assumed that T c = 300 K, 
T h = 400 K, = 0.054 V, and T R = 200 K. 



17 




20 40 j 60 80 100 



